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We measure the number of mechanically stable states of volume fraction of a granular 

assembly under gravity. The granular entropy S{(f>) = log vanishes both at high density, at 
(f) ~ (j>rcp, and a low density, at ~ (l>rvip, where 0rti!p is a new lower bound we call random very 
loose pack, cprip is the volume fraction where the entropy is maximal. These findings allow for a 
clear explanation of compaction experiments, and provide the first first-principle definition of the 
random loose volume fraction. In the context of the statistical mechanics approach to static granular 
materials, states with < (f)rip are characterized by a negative temperature. 



If you pour grains in a container, you will find their 
density to be always enclosed between an upper and an 
lower bound, known as the random close 4>rcp and the 
random loose (prip volume fractions. While (prep has been 
deeply investigated [H , (prip has received considerable less 
attention, and at the present time it has not a precise 
definition, except for the implication that it represents 
the loosest possible, random packing that is mechanically 
stable Here we show that (prip is actually the loos- 
est possible random packing which is mechanically sta- 
ble one can achieve by pouring grains, and that there are 
actually many mechanically stable states of volume frac- 
tion (j) < (prip- We find that the entropy S{(j)) = logr2(0), 
where ^{4>) is the number of disordered mechanically sta- 
ble states of volume fraction <f>, goes to zero at iprcp and 
at (f>rvip < (f'rip-, where (f>rvip is a new lower bound we 
call random very loose volume fraction, (prip turns out 
to be the volume fraction where the entropy is maximal. 
We give an interpretation of both (prcpi 4>rip and (prvip in 
the context of the statistical mechanics of dense granular 
materials. 

To measure the granular entropy, we have considered 
a system of = 20 disks of diameter D under gravity, 
placed in a box with horizontal dimension Ix = 4D, where 
we use periodic boundary conditions. First, we have con- 
sidered a coarse grained lattice model of the system, as 
described in the supplementary materials available on- 
line 0] . The lattice model correspond to a discretization 
of the phase space of the system in cells of size \ as 
we use 51 = D/b, the number of cells is approximately 
[Ix/SVf'^ — 20'"'. Despite this huge number, the entropy 
of the lattice model can be computed exactly, as it is pos- 
sible to reduce the exploration of the phase space to only 
those cells that correspond to mechanically stables state 
which, for the system we are considering, are slightly 
more than 10^. The number of stable configurations at 
each volume fraction 0, ili, ((/>), is shown in Fig. [1^. As 
the lattice states are not actually mechanically stables on 
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FIG. 1: a Number of mechanically stable configurations of 
TV = 20 disks on a lattice, and on the continuous for different 
values of the interparticle friction coefficient, b The volume 
fraction where is maximal approaches the estimated 

value of 4>rip as the friction coefficient increases. 



the continuous, spans unusual small values. In princi- 
ple, the entropy of the system on the continuous could be 
determined investigating that of the lattice model as 51 
decreases; yet, as the number of configurations increases 
as 5l~'^^ , it is clear that this approach cannot be carried 
out. 

Even though the configurations stable on the lattice 
arc not stable on the continuous, they provide a sam- 
pling of the phase space of the mechanically stable states 
of the continuous system. In fact, as schematically de- 
picted in Fig. [2l each configuration stable on the lattice 
is close to energy minima, which can be reached using 
the lattice generated configuration as starting point of 
Molecular Dynamics (MD) of sedimentation under grav- 
ity. We now make an ansatz, which will be justified a 
posteriori, and assume that the lattice generated config- 
urations sample uniformly the phase space of the con- 
tinuous system. Under this assumption, it is possible to 
determine the entropy of the continuous system (up to 
a constant factor) by performing the MD simulations of 
sedimentation, and counting how many minima occurs in 
the volume fraction range (jj, + 5(j). We have performed 
these simulations using a standard model for the grain- 
grain interaction and considering different values of 
the friction coefficient /i. The difference between this ap- 



2 




Coordinates 



FIG. 2: Schematic landscape of tlie energy as a function of 
tlie configurational coordinates. The red disk represents a 
random configuration, used as starting point of energy mini- 
mization procedures to find mechanically stable states (energy 
minima). The blue disks represent the 'almost stable' initial 
configurations we use as starting points of our minimization 
procedure (see text). 

proach to explore the phase space of the mechanically 
stable states and previous ones @,@,0,[1 in only in the 
choice of the initial conditions, as illustrated in Fig. [21 

The number ^^{(j)) of mechanically stable configura- 
tions of volume fraction cj), when the static friction co- 
efficient between the grains is ^, is therefore estimated 
as: 

nM) = I d4>'^L{cb')PM' ^<l>)~ (1) 
- I J deb' J drnL{(i>')nL{cb")PM^'l^" ^<l>)- 

Here the first term accounts for the probability that the 
sedimentation procedure transforms a lattice-stable con- 
figuration of volume fraction 0' in a stable configuration 
of volume fraction </>, the static friction coefficient being 
fi. The second term takes into account the possibility 
that two different lattice-stable configurations actually lie 
in the same energy minimum, and assures that we don't 
count an energy minima more than once. ((/)', 0" —> (j)) 
is the probability that an almost stable configuration 
of volume fraction 0', and one of volume fraction (f>" , 
give rise to the same mechanically stable state of vol- 
ume fraction 0. By performing simulations of sedimen- 
tation using as initial states some of the lattice gener- 
ate configurations, we have estimated Pfi{4>' ~* 'P) &nd 
Pp(0', 4>" 4>) [13. Precisely, we have performed a total 
of 2 10** simulations for each value of /i, each simulation 
ending when the kinetic energy becomes zero within nu- 
merical errors; when this is the case, the system is always 
found in an energy minimum as all of the eigenvalues of 
the Hessian are negative, and therefore no farther relax- 
ation occurs. Figure[T^ shows 0^((/)) for ^ = 0.05, 0.2, 1.0 
and 2.0. For all values of /z, the entropy 5(0) — \og^{(t>) 
has a maximum, and decreases smoothly for smaller and 
larger values of 0. At high 0, Vl{(f)) does not vanishes 



at (f>rcp — 0.82 as we observed the formation of crys- 
talline structures, as usual in 2D dimensional systems. 
The entropy should be an increasing fraction of /i, and 
we observe that this is actually the case but at very high 
density, where the limit of the coarse graining procedure 
used in the lattice model shows up \u\ . 

The value of the volume fraction j^™^'^ where the en- 
tropy is maximal depends on the friction coefficient /i. 
A plot of (j)"^'^^ as a function of yu, as shown in Fig. [T|3, 
reveals that, at high friction, (f)™^^ approaches the limit- 
ing value 0.772. This is a common estimate of (prip for 
frictional systems [l, . Accordingly, we suggest that 
4)rip can be defined as the value (j)^^^ where the entropy 
is maximal. This definition re pro duces the ^ dependence 
of (j)rlp^ previously suggested and does also explain 
compaction experiment as discussed below. These evi- 
dences clarify that there are many disordered mechani- 
cally stable states of volume fraction (f> < <f>rip , and sug- 
gest to introduce a lower bound for the volume fraction 
of a mechanically stable disordered assembly of particles, 
we term random very loose volume fraction, whose esti- 
mation in 2D is 4>rvip — 0.675. To stress that the pres- 
ence of mechanically stable states with very low volume 
fraction is not a finite-size effect, we have also devised 
a procedure to construct low density states with an ar- 
bitrary number of disks, described in the supplementary 
information [3|; a web-applet based on this procedure is 
also available [18j . 

Our results for the entropy are based on the assump- 
tion that the configurations stable on the lattice provide 
a uniform sampling of the minima of the continuous sys- 
tem. While this assumption could be in principle vali- 
dated by investigating lattice models with a small value 
of 51^ we note here that the assumption works, at least 
as a first approximation, because it allows to reproduce 
all expected features of the entropy, as well as to explain 
other experimental evidences found both in 2 dimensions 
and in 3 dimensions. In fact, the shape of the entropy we 
have determined 1) decreases as (j) varies between <f>rip 
and (prep as recently found [H, 2) increases with 
the friction coefficient as expected, 3) explains the de- 
pendence of 4>rip on the friction coefficient, 4) explains 
compaction experiments and why states with volume 
fraction smaller than (prip have not been experimentally 
observed (see below). Moreover, the non-monotonic be- 
havior of the entropy shown in Fig. [T^ is expected from 
the simple consideration that the number of mechanically 
stable configurations of our system can be expressed as 
^{(f) = U,gas{4')^{<j)), where f2gas(0) is the number of 
configurations of volume fraction (p without taking into 
account the stability criterion, and H((/)) is the fraction of 
configurations of volume fraction (f> which are stable. As 
4> decreases, ^gas[4') exp(l/0) exponentially increases 
because there is more space where to put the disks, while 
n(0) decreases. Therefore their product, as well as the 
entropy 5(0), must have a maximum as we have found. 

We now show that the identification of (prip with the 
volume fraction where the entropy is maximal allows to 
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FIG. 3: a The main panel (inset) shows the probabihty P{4>i — > (j>f) (.P[4>f > 4'\4>i)) that a sedimentation procedure which 
starts from an unstable state configuration of volume fraction <j)i ends in a stable state of volume fraction (fjf (greater than (j)). 
b Evolution of the volume fraction in a sequences of taps, c The volume fraction reached at stationary in a tapping dynamics 
as a function of the intensity of the tap 5(j>- d Equation of state in the statistical mechanics approach to granular media. 



explain compaction experiments 0, 9, as well as why 
states with (j) < (j)rip have not been previously observed. 
To this end, wc first consider what is the probability dis- 
tribution P((/)i —f (pf) that a random unstable state of 
volume fraction (pi becomes a mechanically stable state 
of volume fraction 0y in a sedimentation procedure. Ne- 
glecting any dynamical effect, an assumption which is 
reasonable when the sedimentation procedure is very 
slow, we assume that all mechanically stable states with 
volume fraction 4>f < 4>i have the same probability to 
be found via the sedimentation procedure. Accordingly, 
P{4>i — > (pf) is simply proportional to 17(0/): 

while P{4'i = if > (/)/, as during a sedimenta- 

tion process the density increases. 

To calculate this and other quantities we approxi- 
mate in the following the granular entropy with S'(0) = 



\ogri{(p) {(f) — (prvip)" {4>rcp — <PY (s^e dashed line in 
Fig. [1^), but different approximations give rise to very 
similar results. We show in Fig. [3^, for different values 
of both the probability distribution P{<f)i Pf) and 
the probability distribution that p>f > (p, P{4>f > 4>\ Pi) 
(inset). If (pi < (prip, P{4>i ~* 4>f) has a maximum at 
(pf — (prip- If Pi > 4>rip, on the contrary, the maximum 
of P{pi (f>f) is at (f>f ~ (pi. These results indicate 
that in usual sedimentation procedures packs of volume 
fraction smaller that (j)rip are not observed just because 
of entropic reasons: the number of states with (j) < (jirip 
is exponentially too small with respect to the number of 
packs of volume fraction (p ~ (prip ■ 

Eq. [2] could be used to describe compaction experi- 
ments, where a granular system is subject to a sequence 
of taps. In a single tap, a mechanically stable state of vol- 
ume fraction 0„ first becomes an unstable state of volume 
fraction </>„ — Sp, due to the mechanical energy which is 
given to the system, and then settles down in a new me- 



chanically stable state of volume fraction (j>n+i- 6(j) > is 
related to intensity of the tap: the stronger the tap, the 
larger S(j). After many taps it is observed that the volume 
fraction reaches a steady state 4>{64>) which depends on 
the intensity of the tap Fig.[3)D shows that this 

kind of behavior is reproduced by considering (/)„ — 5(f) ^ 
4>i and 0„+i ((/)„, (5(/)) = / (j)fPii4>n - Scf) 4>f)d(j)f. As 
the number of taps increases the volume fraction reaches 
a steady value 4>{S(t)), which is the fixed point of the re- 
currence relation 4>n+i ~ 0„+i(0ri, The dependence 
of (j){S(j)) on 5(j), which is shown in Fig. [St, reproduces 
the dependence of the volume fraction reached at sta- 
tionary on the intensity of vibration in compaction ex- 
periments d, H, 0, [1| ■ The identification of (j)rip with the 
volume fraction where the entropy is maximal is there- 
fore consistent with previous works which have consid- 
ered (firip as being the loosest pack achieved via sedimen- 
tation or tapping procedures [3, 0, [HI • 

The results presented here have an interesting inter- 
pretation in the context of the statistical mechanics ap- 
proach to static granular media In this theoret- 



ical framework, one considers a granular system sub- 
ject to a dynamics allowing for the exploration of the 
mechanically stable states, and assumes that all states 
of given volume V are visited with the same probabil- 
ity. For hard spheres under gravity the energy E is 
dominated by the gravitational energy, and therefore 
V (X E (X l/(/>. This hypothesis has been proved numeri- 
cally at ~ (f>rcp for a system under shearing [isj , and in 
the density range (jirip < 4> < 4>rcp for particles driven by 
flow pulses Accordingly, one can introduce a granu- 
lar temperature, known as configurational temperature 
or compactivity, defined as T'J^j: = (iconf = dS/dE. 
The dependence of (j) on /?con/ (see Q for calculation 
details) is shown in Fig. [SJi, and clarifies that in this the- 
oretical framework 0rcp, i'rip and 4>rvip are defined as fol- 
lows: cj)rcp = (t){l3conf +Oo), 0r-ip = <t)iPconf = 0), and 

4>rvip = (Pif^conf — oo). The density range that is not 
usually observed, with volume fraction enclosed between 
(prvip and (firip, is that where the granular temperature is 
negative. 
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Supplementary information, 
Methods 

Hard disks on a lattice 

In order to uniformly sample the phase space of the mechanically stable configurations of an hard disk sys- 
tem under gravity, we have used as starting point of a sedimentation procedure configurations which are 'almost 
stable', meaning that they are close to an energy minima. We have generated these configurations using a lattice 
model obtained by coarse graining an hard disk system, the 'almost stable' configurations being defined as those that 
are stable according to the stability rules defined on the lattice. 

In the model, the presence of a grain in a given lattice site hinders the occupancy of nearby lattice sites due 
to steric constraints. For instance, the grain shown in Fig. |4] forbids the occupancy of the lattice sites marked with 
a star (as well as of all the lattice sites covered by the disk). As we work under gravity, we consider a grain to be 
stable if it is supported by both a grain on its left and by a grain on its right, i.e. if at least one of the positions 
marked with L and one of the positions marked with R in Fig. |4] are occupied. We also consider the possibility that 
two grains A and B cooperate to their stability. Precisely, if grain A is in position and grain B in position 

{i + D,j), and both one of the i-positions of grain A and one of the i?-positions of grain B are occupied, then both 
grains are mechanically stable. Finally, a grain resting on the bottom or on the first layer (j = 0, 1) of the lattice 
is also considered as mechanically stable. A configuration is mechanically stable if all of its grains are mechanically 
stable. The stability rules implies that, if a grain at a given heigh h is not mechanically stable, then it is not possible 
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FIG. 4: Disks on a lattice. If a disk is present in the system as in figure, then all of the lattice sites covered by the disk as well 
as the lattice sites marked with * cannot be occupied by other disks due to steric constraints. The positions marked with L 
and R are used to define the mechanically stable configurations in the text. 

to make it stable by adding grains at an height h' > h. By exploiting this feature it possible to devise a procedure 
to selectively generate all of the mechanically stable configurations of this model, without the need to generate the 
unstable ones. As we want to evenly sample the phase space, we have considered a system of = 20 disks on a 
square lattice with horizontal size AD and lattice constant D/5, for which it is possible to generate all of the stable 
configurations. A typical lattice generated configuration, together with the configuration reached after sedimentation, 
is shown in Fig. (5] 

It is not possible to investigate significatively larger values of A", as the number of mechanically stable configurations 
grows exponentially with A^, as shown in Fig. [6l 
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FIG. 6: The number of lattice-stable configurations increases exponentially with the number of grains. 



Supplementary information, 
Methods 

Generation of very loose mechanically stable random states 

Via the combined use of tiie lattice model to generate 'almost stable' configurations, and of Molecular Dy- 
namics simulations of sedimentation under gravity, it is possible to uniformly sample the phase space of an hard 
disks system, and to generate mechanically stable states of with (f> < 4>rip- This approach has the advance of allowing 
an uniform exploration of the phase space, but is limited to systems with a small number of particles. 

Here we present a procedure to generate mechanically stable states with a very small value of (f> with an arbitrary 
number of particles, in the limit of hard particles with very high friction. A web-applet based on this procedure is 
availabe (http://smcs.na.infn.it/rvlp/rvlp.html). Figure [7] shows a pack generated with this procedure, with volume 
fraction (calculated in the area delimited by the two horizontal lines) 4> — 0.746. 




FIG. 7: A mechanically stable loose with volume fraction ~ 0.746 < (jirip — 0.772. 

The procedure is based on a growing mechanism, where disks of diameter D are added on top of previous placed 
disks. Given a mechanically stable state, the number of particles can be increased by adding n disks at the same time. 
If n = 1, then there is a finite number of positions where the new grain can be placed. These positions can be easily 
found from geometrical considerations and are show as blue shaded disks in Fig. [T] For n > 1, there is an infinite 
number of mechanically stable configurations the disks n can be arranged. We have considered, for simplicity sake, 
only n = 2, and we have placed the two grains according to the following steps: 

1. Choose a disk A on the bed surface and an angle 6, —it < 6 < tt. 

2. Place a new disk a in contact with disk A, the line joining their center making an angle 9 with the positive 
vertical direction. If the grain overlaps other grains, eliminate the new grain and go back to 1. 

3. Determine the positions where it is possible to add a single grain, the grain being in contact with grain A. If 
there are no positions, go back to 1. 

4. Place a new grain (3 in one of the previously determined positions. The grain will touch a grain B. If grain a 
and P are not stable, eliminate them and go back to 1. 

In the limit of hard disks and infinite friction we are considering, the stability of grains a and /3 can be easily 
determined form geometrical considerations, as illustrated in Fig. [S] For instace, grain a is stable if xa < Xa < xb 
and yp > yp, where P is the point where the line joining grains A and (3 and the vertical line passing trough the 
center of grain a intersect. A similar condition holds for grain /3. The resulting bridge is stable if both grain a and 
grain /3 are stable. 




FIG. 8: Stability of a bridge. 



The hypothesys of infinite friction can be relaxed after the packs has been build. In fact, if the pack is used as starting 
point of soft core Molecular Dynamics simulations of relaxation under gravity, one observes a small compactification, 
related to the softness of the disks, during which tangential forces build up. Once the pack becomes stable, one can 
determine the the minimum friction coefficient required for the pack to be stable, by measuring the normal Fn^ and 
the tangential f/*-* forces at each contact i. The minimum friction coefficient is: /i = max^ |. 



Supplementary information, Methods 

The equation of state 

The equation of state shown in Fig.4d, expressing the dependence of </> on the inverse temperature is easily 
determined as 

We use for the entropy the expression 

which describes the numerical data (see Fig. la) with A = 7, 4>rvip = 0.675, 4)rcp = 0.87, a = 0.28, (3 ~ 0.15. The 
energy E is given by the gravitational energy of the system, 

E = ^Nmgh, 

where = 20 is the number of disks, m the mass of each disk, g the gravitational field, and h the height of the 
system, h is related to the volume fraction cf), as 

, Nv fD^ ^ 



where D is the diameter of the disks, and Ix = 4_D the horizontal size of our system. 
From the above expressions, we obtain: 



dS , 



and 



